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ABSTRACT 

At the Rocky Flats Wind Systems Center, several 
different tower dynamics analysis methods and 
computer codes are used to determine the natural 
frequencies and mode shapes of both guyed and 
freestanding wind turbine towers. In this paper 
these analysis methods are described and the 
results for two types of towers: a guyed tower 

and a freestanding tower are shown. The advan- 
tages and disadvantages in the use of and the 
accuracy of each method are also described. 

INTRODUCTION 

The accurate prediction of tower vibration fre- 
quencies and mode shapes Is important 1n avoiding 
unwanted vibration problems. At present, there is 
a variety of structural dynamic analyses covering 
a range of complexity and application. In this 
paper, some of the existing analyses and corres- 
ponding computer codes will be examined in order 
to determine those which can be of use to the 
SWECS Industry. 

A tower dynamics supporting research and technol- 
ogy project has been conducted at the Rocky Flats 
Small Wind Systems Center. The objective of this 
project has been to determine those analyses which 
are simple to use but give adequate results com- 
pared to test results. This paper will present 
some of th^ simpler tower dynamic analyses, their 
correct use and accuracy. The simpler analysis 
methods will be presented in order of increasing 
complexity and accuracy. The theory and use of 
the iiethods will first be described, as well as 
the accuracy of results for towers with various 
mass and stiffness distributions. The results 
from these analyses will then be compared to the 
test results for two types of towers at Rocky 
Flats: 1) a guyed tower with uniform mass and 

stiffness, and 2) a freestanding tower with 
uniform mass but tapered stiffness. It will be 
shown that accurate determination of the bending 
frequencies for towers with tapered stiffness is 
mre difficult, using the simpler analyses. This 
will require the designer to use^ a more complex 
analysis, such as SAPIV. It should be emphasized 
that this paper deals specifically in determina- 
tion of tower bending frequencies. The case of 
torstona! frequencies or coupled bending torsional 

*Now”with Energy Sciences, Inc., Boulder Colorado. 


frequencies have not been analyzed, using the 
methods of this paper. The results for these 
frequencies may be the subject of a later paper. 

THEORY AND USE OF THE METHODS 
The Rayleigh Quotient 

The Rayleigh Quotient forms the basis for some 
approximate methods: the Rayleigh Method and 

the Rayleigh Ritz Method, both to be described 
here. The system to be analyzed is shown in 
Figure 1. The freestanding tower has variable 
mass and stiffness distributions pA(x) and 
EI(x), respectively. A concentrated mass, M, is 
attached at the tower top at a height, L, above 
the ground. 



Figure 1: System to be Analyzed. 

By assuming the deflection v(x,t) in the form: 
v(x,t) <j»(x) sinpt 

and equating the tower potential energy at 
maximum displacement to the kinetic energy at 
minimum displacement, the Rayleigh Quotient can 
be expressed: 
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L ? 

/ EKx) [♦"(x)]'= dx 

p2(+) = I 

/ pA(x) [Mx)]^ dx + M[^(L)]^ 


From Table I it can be seen that the accuracy of 
results from Equation 1 diminishes for towers 
with higher stiffness taper rates and larger tip 
masses because the assumed first mode shape 
(|)(x), used in the above analysis, approximates 
the exact mode shape less accurately. 


When <fi Is an exact vibration mode of the tower, p 
has the value of the corresponding exact natural 
frequency; also, p^ is stationary with respect to 
variations in <|) at each of these points {Ref. 3). 
This expression can be expected to give a good 
approximation to the frequency if a good approxi- 
mation to the mode shape is used to evaluate 
Equation 1. A simple analysis can be formed for 
the special case of a tower with uniform or 
linearly tapered mass and stiffness distributions. 
In order to find the first bending frequency. 

A Simple Analysis 

For a freestanding tower with linear mass and 
stiffness distributions of the form; 

EI(x) = EIq (1-0 x/L) 

pA(x) = pAq (1 - Y x/D* 

the approximate shape: 

4 3 2 

<^'(x) = (f) - 4 (^) + 6(^) ** 


can be used to evaluate Equation 1 with the 
result: 


El ^ 28.800 - 4.8003 

pAqL^ (2.311 - 1.854t) + 9.00p 


( 2 ) 


TABLE I: COMPARISON OF RESULTS OBTAINED FROM 

SIMPLE METHOD (W.) TO EXACT RESULTS 
(Wp) FOR THE FIRST BENDING FREQUENCY 
RATIOS, FOR VARIOUS TAPERS AND TIP 
MASS RATIOS, 

Uniform Stiffness, Uniform Mass: 0=0., y = 0 



w 

"e 

«s 

% Error 

1.0 

1.56 

1.60 

2.6 

2.0 

1.16 

1.19 

2.6 


Uniform Stiffness, Tapered Mass: 0=0., y = .9 


y 

W 

"e 

«s 

% Error 

1,0 

1.68 

1.73 

3.0 

2,0 

1.21 

1.24 

2.5 


Tapered Stiffness, Uniform Mass: 0 = .9, y = 0 


y 

We 


% Error 

1,0 

1.33 

1.47 

10.5 

2.0 

1 0.99 

1.10 

11.1 


Limitations of the Method 


where m is the ratio M/pAqL. 

In those cases where exact frequencies for a 
cantilever beam with a tip mass have been calcu- 
lated, the results from Equation 2 can be com- 
pared directly to give an indication of the 
accuracy. Table 1 shows the dimensionless 
frequency ratios: 



This simple method is applicable to freestanding 
towers for mass and stiffness distributions 
varying approximately linearly. Only the first 
bending frequency can be found from Equation 2, 
although the second bending frequency might be 
estimated if the analyzer could evaluate Equa- 
tion 1 with a good second mode shape approxima- 
tion. Equation 2 is useful for obtaining a 
rough estimate of the first bending frequency of 
freestanding towers, without the use of a 
computer. 

A Rayleigh Computer Code 


obtained from Equation 1 compared to exact 
results calculated, using a power series method 
(Ref. 8). 


*EIo and pAq are the stiffness and mass 
distributions of the tower at the base. The 
case 0 = 0, Y = 0 is that of a uniform tower. 

**The assumed mode shape should satisfy some or 
all of the cantilever tower boundary condi- 
tions. Here: 4>(o)= 0 and <|)'(o) = 0 are satis- 

fied. The assumed mode shape should at least 
satisfy the base boundary conditions. 


A program developed at MIT, Program Rayleigh 
(Ref. 4), uses the Rayleigh procedure to evalu- 
ate the first bending mode for freestanding 
towers with nonlinear mass and stiffness distri- 
butions. The tower can be divided into N equal 
segments and constant values of mass and stiff- 
ness are input along each segment, Simpson’s 
integration is then used to evaluate the inte- 
grals in the Rayleigh Quotient. 

Both the Simple Method and the Rayleigh Program 
have the disadvantage that the user must input 
an assumed mode shape. The accuracy of the 
corresponding results depends on this one shape. 
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In a ipethpd to be described next, a linear com- 
bination of assumed shape functions is used to 
form the function ^(x) for use in the Rayleigh 
Quotient. This method is particularly useful for 
guyed towers, because the mode shapes are harder 
to approximate with a simple function, such as 
that used in the simple method above. 

THE RAYLEIGH RITZ PROCEDURE 


Jhe Rayl eigh Ritz method Involves using a set of 
assumW shape functions i>-j(x) and combining 
them to form the nx)de shape: 

+ A2i|^(x) + . , . + (3) 


The Aj's are constants to be determined and the 
^ s are a set of linearly independent func- 
tions, each satisfying some or all of the boun- 
dary conditions of the tower.'*' 

Per a guyed tower, Rayleigh's Quotient can be 
expressed!^ 


J EKx) [$"(x)]^ dx + Kc[ 4 .(a)]^ 

* =t 

- / pA(x) [♦(x)]2 dx + H[.j.(L)]^ 

0 


(4) 


where Kc is the guy wire stiffness coefficient 
and a is the guy wire attachment height, shown in 
Figure 2. For three guy wires spaced 120*^ apart, 
Kc can be shown to be (Ref. 6); 


Kc - ^ ^ — ^cos^B. (5) 


While for four guy wires spaced 90® apart, the 
factor 3/2 in Equation 5 is replaced by a 2, 

These values of Kc are valid when the guy wires 
have been tensioned sufficiently so that there is 
no coupling between the guy wires and the tower 
(Ref. 1). 


Substitution of the mode shape form in Equation 3 
into Equation 4 gives: 


El 



Where: 


5 1 N *j ‘u 


J 1 


(6) 


1 



II 

o 


"( Ode + Kc ( ea) <i>j( ea) 


1 




= / h(d 
0 


Ode + ( 1) >i>j { 1 ) 


i, J = 

1..... N 


and 

C, fU), 

h(5), 5a 

, Rc and \i are the dimen- 


slonless distance, stiffness, mass distribution, 
guy level attachment, guy stiffness coefficient, 
and tip mass ratio, respectively, as defined in 
the nomenclature. 



Figure 2: Guyed Tower to be Analyzed. 

It is known that the Rayleijh^ Quotient is sta- 
tionary with respect to” WrTafti o n s i n <f»(x), when 
(j)(x) is an exact mode shape for the tower. This 
can be shown to require that: 

Det [d,-j - W^Tj] = 0 (Ref. 3). 

The problem is thus reduced to finding the H 
values of which make the determinant of this 
N X N matrix equal to zero. The N values of W 
give estimates for the first N frequencies of 
the tower. 


Simple Analy si s 

A simple two-mode analysTs for freestanding 
towers can be performed by using two approxi- 
mating functions of tbe~Tbrm: 


111(5) 

= 5“ - 

4C^ + 


i(2( 5) 

= 35= - 

- lO^** + 10 


The integrals d-jj 

and b^-j then 

are: 

dll = 

28.800 

- 4.800ft 


di2 = 

24.000 

- 6.857ft 


d22 = 

34.286 

- 12.857ft 


bii = 

2.3111 

- 1.8540T + 9 

.00^4 

bi2 = 

2.0698 

- 1.6997Y + 9 

.00^ 

b22 = 

1.8817 

- 1.5732Y + 9 

.00]i 


when det [d^j - is evaluated, the 

polynomial glW^)^ - qjW^ + i= 0 must be solved 
where: 

9 = biib22 - 

q = d]^^b22 ^11^22 " ^^12^12 

I = diid22 - 


*Thc function^ Should satisfy at least the 

geometric boundary conditions (Ref. 3). 
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Table II shows the results for various values of 
tip mass ratio and taper rates. This shows that 
the simple two-mode method becomes inaccurate for 
the second bending frequency ratios of towers 
with large stiffness tapers. This method gives 
more accurate results for the first bending fre- 
quencies compared to the results of the previous 
sections. 

TABLE II: EXACT (W^e^ APPROXIMATE (W^a) 

FIRST AND SECOND BENDING FREQUENCY 
RATIOS FOR A FREESTANDING TOWER, 
TWO-MODE RAYLEIGH RITZ PROCEDURE, FOR 
VARIOUS TAPERS AND TIP MASS RATIOS. 

Uniform Stiffness, Uniform Mass: 6=0, r 0 



le 

“la 

% Error 

“2e 

“2a 

% Error 

1.0 

1.56 

1.56 

0.0 

16.25 

17.61 

8.4 

^0 

1.16 

1.16 

0.0 

15.86 

17.29 

9.0 


Unform Stiffness, Tapered Mass: 3=0, y = ,9 


tl(«) = 5" - + 6?2 

42(5) = - IOC- + lOc^ 

43(5) = 25® - 6^5 + 5?- 

44(C) = c'® - 2.5c® + 1.6075®. 

The program can be used for linear or nonlinear 
mass and stiffness distributions by dividing the 
tower into N sections (not necessarily equal) 
and inputting constant values of stiffness and 
mass along each section. Numerical integration 
techniques are then used to evaluate the Inte- 
grals in the Rayleigh Quotient. Table III shows 
the improvement for the first and second bending 
frequency results, compared to the simple two- 
mode method of Table II. 

TABLE III: EXACT (W^g) AND APPROXIMATE (W^^) 
FIRST AND SECOND BENDING FREQUENCY 
RATIOS FOR A FREESTANDING TOWER. 
FOUR-MODE RAYLEIGH RITZ PROCEDURE FOR 
VARIOUS TAPERS AND TIP MASS RATIOS. 


p 

le 

“la^ 

% Error 

“9 

2e 

“2a 

% Error 

1.0 

1.68 

1,70 

0.6 

22.37 

23.56 

5.3 

2.0 

Taperec 

p 

1.21 
i Sti1 

le 

1.21 

Ffnes! 

w. 

1 a 

0.0 

5, Uni for 
% Error 

22.17 
m Mass: 

2e 

23.00 
6 = 

3.7 

.9, Y = 0 
% Error 

1.0 

1.33 

1.37 

3.0 

11.81 

15.36 

30.1 

2.0 

0.98 

1.02 

4.1 

11.53 

15,06 

31.0 


This simple hand analysis can also be used for 
guyed towers, using these two approximating func- 
tions and including the guy wire stiffness 
effects in the terms d^j. As an example, a 
single guyed tower with the dimensionless stiff- 
ness ^c = 200 and guy level ^ = .8 was examined 
for the first two bending frequencies. It was 
found that W^ « 20.16 and W9 ~ 23.43, These 
same results were calculated with a computer 
program, to be described in the next section, 
using a four-mode Rayleigh Ritz procedure with 
the results: W^ ~ 19.82 and W2 « 22,45. 

These results show that the simple two-mode 
Rayleigh Ritz procedure gives fair results for a 
guyed tower al so . 

As will now be shown, the use of more approxima- 
ting functions will improve the results of 
both the first and second bending mode frequen- 
cies and will also give results for the higher 
mode frequencies. 

The Rayleigh Ritz Computer Program 

The Rayleigh Ritz method can be used to determine 
any number of modal frequencies; however, the 
computational complexity increases greatly after 
the first two or three modes. For this reason, a 
computer code has been developed at Rocky Flats 
using four approximating functions of the form: 


Uniform Stiffness, Uniform Mass: p = 0, y = 0 


p 

le 

1 a 

% Error 

j 

Wo 

2e 

Wo 

2a J 

% Error 

1.0 

1.56 

1.56 

0.0 

J 

16.25 

16.26 

0.1 


1.16 

1.16 

0.0 


15.86 

15.88 

0.1 

Unform 

Stiffness, 

, Tapered Mass: 

p = ( 

3, Y = .' 

i. 

^le 

W. 

la 

% Error 


J 

“2a 

% Error 

1.0 

1.68 

1.70 

0.6 


22.37 

22.49 

0.5 

2.0 

1.21 

1.21 

0,0 


22.17 

22.30 

0.6 

Tapered Stiffness, Uniform Mass: 

: P = 

•9, r = 

p 

“le 

la 

% Error 


“2e 

“ 2 a 

1% Error 

1.0 

1.33 

1.34 

1.0 


11.81 

12.75 

1 8.0 

2,0 

0.98 

.93 

1.0 


11.53 

12.45 

8.0 


The program can also be run for guyed towers, 
when the guy wire tension is large enough so 
that no coupling between the guy wire first mode 
and tower first mode occurs (Ref. 1). 

Limitations of the Rayleigh Ritz Method 

Both the simple hand method and the four-mode 
Rayleigh Ritz program can be used to calculate 
first and second mode bending frequencies for 
freestanding and guyed towers. The simple 
method gives poor results for the second bending 
frequencies of towers with high stiffness 
tapers. The more complex four-mode method can 
be used for towers having more complex nonlinear 
mass and stiffness distributions, and gives 
better results than the two-mode method, for 
highly tapered towers. The disadvantage of this 
method is that a computer must be used. 
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A code for a TI-59 programmable calculator has 
been developed to handle towers with nonlinear 
mass and stiffness distributions* 

Program TUSF (Turbine ^stem £requencies) 

The two-nK)de Rayleigh Ritz procedure covered 
above can only be used for towers which have 
approximately linear mass and stiffness distribu- 
tions. For towers with nonlinear distributions, 
this method v^uTd be hard to use, because the 
integrals in the Rayleigh Ritz Method would be 
too hard to evaluate. 

To circumvent this difficulty, a method similar 
to the Rayleigh Ritz procedure has been developed 
for towers with nonlinear mass and stiffness 
(Ref. 6). The tower can be divided into H sec- 
tions and constant values of mass and stiffness 
are input along each sigment. Numerical integra- 
tion techniques are then used to find the natural 
frequencies (Ref. 7). 

In th e techniq ues described previously, the 
nacelle and rotor were modeled as a single lumped 
ma^ ^ Uie tower top. This program mludes the 
effects of rotor moments of inertia and nacelle- 
rotor C.G. location on the system frequencies. 

In addition, the effects of rotor spin rate on 
the natural frequencies are taken into account. 

The program treats the tower as a flexible member 
and the nacelle and rotor as rigid bodies. It 
will be shown in the section on conparisons of 
analytical results to test results that the 
accuracy of results from this program are similar 
to the accuracy of results from the Rayleigh Ritz 
procedure, I.e., the results for the second 
bending frequencies are less accurate than the 
first bending frequencies. 

In the approximate methods discussed above, an 
Infinite degree of freedom structure is modeled 
as a beam having finite degrees of freedom. It 
can be shown that this causes the frequency 
estimates from these methods to be higher than 
the exact values (Ref. 6). A method using more 
degrees of freedom (such as SAP IV) gives more 
accurate frequency estimates. 

In the next section, the test procedures and 
results for two types of towers will be present- 
ed. In the section on comparisons of analytical 
results to test results these simple methods, for 
these towers, will be compared. The results from 
more accurate structural codes (such as SAP IV) 
will also be shown, A recommendation as to which 
tower types can be analyzed, using the simple 
method, will also be given. 

TEST PROCEDURES AND RESULTS 
Objectives 

This section deals specifically with the test 
methods and results for a Rohn 25G single guyed 
tower and a Rohn SSV freestanding tower. Impact 
tests were performed on each tower to determine 
the predominant modes of vibration and to compare 
the test results with various analytical methods. 


Test Description and Theory 

A technique currently being used at the Rocky 
Flats Test Center, for dynamic testing, is 
commonly known as impact testing. The structure 
can be excited with an impulse. This can be 
accomplished with the use of a hammer with a 
load cell attached, thereby exciting the struc- 
ture with a known input, as shown in Figure 3. 



Figure 3. 

With the load cell, the input force can be 
accurately measured; the response of the struc- 
ture can be measured with the use of an acceler- 
ometer. Provided the input and response signals 
are fourier transformable (Ref. 2), the frequen- 
cy response function can be computed. 

In practice, better results are obtained by 
computing frequency response functions as a 
ratio of the cross spectrum between the Input 
and output, to the power spectrum of the input 
(Ref. 2), This 1s useful In eliminating the 
effects of noise on the input and output signal 
measurements. If the input and output noise is 
noncoherent, the effects on the cross spectrum, 
involving the noise, will yield zero. 

Test Specimens and Test Results 
A Rohn 25G Single Guyed Tower 

Figure 4 shows a Rohn 25G 40-ft guyed tower. 

The tower has constant mass and stiffness 
distributions as given in Table IV. The tower 
stiffness distribution was determined by first 
finding the moment of inertia of a tower cross 
section about the centroid, as shown in Fig- 
ure 5. Because the distance between the tower 
legs remains constant, for various heights, the 
moment of inertia and thus the stiffness distri- 
bution remains constant. The mass distribution 
also is constant, since the weight and length of 
each tower section is the same. The effects of 
the cross braces have been neglected in the” 
stiffness determination, but have been included 
in the mass distribution. 


t 
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Figure 4. 

TABLE IV: ROHN 25G PARAMETERS 

Tower Height (L) 40.0 ft 

Stiffness Distribution (El) 3.14 x 10^ Ib-ft^ 


Mass Distribution (pA) 

Tower Top Mass (M) 

Guy Wire Attachment Level 

Guy Wire Cross Sectional 
Area (AM 

Guy Wire Elastic Modulus 

Guy Wire Length { Jt) 

Guy Wire Stiffness 
Coefficient (Kc) 


0.124 Ib-sVft^ 
10.87 Ib-sVft^ 
(a) 35.0 ft 

0.00015 ft^ 

(EM 2.88 X 10^ 1b/ft^ 
48.45 ft 

6365.76 Ib/ft 



Figure 5. 


The Rohn 25G was tested without a machine on 
top, but with the addition of 380 lb on the 
tower top. Table V shows the test results, 
using the impact testing methods. In this case, 
the three guy wires were tensioned sufficiently 
so that each had nearly the same fundamental 
frequency. As can be seen from Table V, the guy 
fundamental frequency was well separated from 
the tower first bending frequency, so that 
resonance between the guy wires and tower was 
not a problem (Ref, 1). 

As will be shown in the next section, the simple 
methods can be used to get good estimates for 
the first bending frequencies for this tower, 

TABLE V: ROHN 25G TEST RESULTS 


First Mode Bending 2,6 Hz 
First Mode Torsional 5.7 Hz 
Second Mode Bending 7.8 Hz 
Guy Wire Fundamental 6,9 Hz 


A Rohn SSV Freestanding Tower 

Figure 6 shows a 60--ft Rohn SSV freestanding 
tower. Unlike the Rohn 25G, the distance 
between tower legs decreases with Increasing 
tower height. The stiffness distribution thus 
decreases from base to height. Figure 6 shows a 
plot of the mass and stiffness distribution for 
this tower. The mass distribution is nearly 
uniform, but the stiffness distribution tapers 
to 0.023 of the value at the tower base. 

The Rohn SSV was tested without a machine on top 
but with a top plate of approximately 30 lb. As 
will be shown In the next section, accurate 
estimation of the first and second bending 
frequencies for this tower Is difficult, using 
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the simple n^thods, because of the large stiff- 
ness taper. Also, the effects of the cross- 
meiTt)ers were neglected In the moment of Inertia 
calculations, which poses further difficulties. 
Table VI shows the test results for the first two 
bending frequencies. 

TABLE VI: ROHN SSV TE_ST RE.SULTS 

First i^de Bending 3.3 Hz 

Second ^tode Bending 12,7 Hz 

COMPARISONS OF ANALYTICAL RESULTS TO TEST RESULTS 

In_thl 5 section the results from the simple anal- 
ysis methods and results from complex methods 
such as SAPIV (Ref. 5) will be compared to test 
results for the two towers presented in the 
previous section. A recommendation of the use 
and accuracy of the simple methods will be made. 

Rohn 25(^ Results Suimnary 

Table VII shows the results for the Rohn 25G. 
ThflimpTe two -mode Rayleigh R1tz method, a hand 
method, gives fatr results for the first frequen- 
cy (pi^), but poor results for the second fre- 
quency (p2)- 1^he utility of this method is 
that it can be used to obtain a rough estimate, 
by hend”c^cuT a ti on for the 71 r s rtendi h g 
frequency of guyed or freestanding towers. This 
method gives less accurate results for towers 
with large stiffness tapers, however. 


" TABLE VI I: ROHN 25G RESULTS SUMMARY 


- Method 

Pi 

Error 

P2 

Error 

Test 

2.6 

0.0% 

7.8 

0% 

R-R-2 Mode 

2.8 

7.7% 

9.0 

16% 

Program TUSF 

2.6 

0.0% 

11.2 

44% 

R-R-4 Mode 

2.6 

0.0% 

8.6 

10% 

MITGUY 

2.6 

0.0% 

8.6 

10% 


Program TUSF, a hand calculator code, gives much 
mre accurate results for the first mode. This 
code can also be used to determine the first 
bending frequency of towers with nonlinear mass 
and stiffness distributions, without use of a 
large computer. 

The Rayleigh Ritz four -mode procedure, a small 
computer code developed at Rocky Flats, gives 
excellent results for the first mode and fair 
results for the second mode, while Program MITGUY 
(Ref. 4), a program utilizing matrix iteration 
methods, gives abovit the same results (to one 
decimal place). 


The 10 percent inaccuracy In the second mode 
results Is thought to be caused because of 
neglect of the cross-members In determination of 
the tower bending stiffness. This difficulty 
will be emphasized In the comparison of results 
for the Rohn SSV, to be shown next, In which use 
of the complex code SAPIV, which models the tower 
cross-ment>ers. Is necessary. 


Rohn SSV Results Summary 

Table VIII shows the results for the Rohn SSV 
tower. The analysis of the Rohn SSV with the 
simple hand methods of the Rayleigh Method and 
Two-Mode Rayleigh Ritz Method Is difficult, 
because the stiffness and mass distributions are 

nonlin ear. For th ese analyses, the m ass 

distribution was as sumed constant, with a value 
equal to the base value. The stiffness 
distribution was approximated as linear, with a 
value of the taper rate 3 equal to 0.98. The 
resulting frequency estimates are too high, 
because the actual stiffness of the tower is 
less than what is given by this straight line 
approximation. 


TABLE VIII: ROHN SSV RESULTS SU?^ARY 


Method 

Pi 

Error 

P2 

Error 

Test 

3.3 

0% 

12.7 

0% 

Rayleigh 

3.8 

19% 

— 

— 

R-R-2 Mode 

3.6 

12% 

22.3 

76% 

Program TUSF 

3.5 

9% 

16.3 

28% 

R-R-4 Mode 

3.4 

6% 

15.5 

22% 

MITGUY 

3.4 

6% 

14.8 

17% 

SAPIV 

3.1 

3% 

13.0 

2.4% 


For towers with nonlinear mass and stiffness 
distributions, the program TUSF and the four- 
mode Rayleigh Ritz procedure can be used to find 
more accurate values of the first bending fre- 
quencies, The tower can be divided Into N 
sections and constant values of mass and stiff- 
ness input along each segment. It was found 
that about 20 segments were needed In order to 
model the large stiffness taper correctly. 

From Table VIII It can be seen that the second 
bending frequency results from: 1) the simple 

methods, and 2) the program MITGUY are 1n error 
by more than 17%. For this reason, the complex 
code SAPIV (Ref. 5) was used to model this 
tower, taking Into account the effects of each 
cross-merrtier. The improvement in the second 
bending frequency Is very evident. 

In the simple methods, as well as MITGUY, this 
tower was modeled as a beam, neglecting the 
effects of the cross-members In the bending 
stiffness. The cross-section moment of inertia 
was calculated, by taking Into account the three 
legs only, as shown In Figure £. The actual 
stiffness distribution of this tower Is more 
complicated than this approximation because of 
the cross-members. 

In the simple methods section It was shown that 
they give less accurate results for a beam with 
large stiffness tapers. For this truss-type of 
tower, this error Is Increased because the 
actual tower stiffness 1s more complicated than 
that of a beam 1n which the cross-members have 
been neglected. 
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CONCLUSION AND RECOMMENDATIONS 

The simple methods can be used to determine good 
estimates of the first bending frequencies and 

rough estimates _of the second f^r’equencie^ of 

towers with uniform mass and stiffness distribu- 
tions, Guyed towers can also be analyzed, using 
these methods. If the guy wire fundamental fre- 
quency is well separated from the tower first 
bending frequency so that guy wire-tower inter- 
action will not occur. 

For towers with high stiffness tapers, or stiff- 
ness d i s tr i buttons which are h ar d to de termine 
accurately, more complex codes such as SAP IV may 
be needed, especially for accurate determination 
of the second and higher bending frequency. 
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NOMENCLATURE 

Guy wire attachment level. 

Guy cable cross-sectional area. 

Guy cable elastic modulus. 

Tower stiffness at station x. 

Tower stiffness at the base. 


a 

A’ 

E' 

EKx) 

Eln 


f(C) 

WiK) 

Kc 

ITc 


L 

i : 

M : 

P • 

t : 

v(x,t): 
W : 

Ai : 


^1 


Y : 
K : 

6 : 

W : 
pA(x): 

pAq : 
♦(x) : 


EI(L5)/EIo 
cA(L 0/ pAp 

Guy wire stiffness coefficient. 

Dimensionless form of Kc: 

^ 3 

Kc = KcL 

nr 

Tower height. 

Guy cable length. 

Lumped mass at tower top. 

Bending frequency (rad/sec) or 
(Hz). 

Time (sec.). 

Tower deflection. 

Dimensionless frequency. 

Coefficients in approximate 
mode shape. 

Functions for approximating 
mode shapes. 

Rate of decrease of linear 
stiffness. 

Rate of decrease of linear mass 
di stribution. 

X/L. 

a/L. 

Angle guy cable makes with 
ground. 

Tip mass ratio: M/pAqL. 

Mass distribution of tower at 
section x. 

Mass distribution of tower at 
tower base. 

Tower mode shape. 
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QUESTIONS AND ANSWERS 


Frum 

Q: 

A: 

Fronts 

Q: 

A; 


A.D. Wright 


: G* Beaulieu 

1) Don’t you believe that direct solution of 4th order differential equation 
would give more accurate mode shape and frequencies? 

2) Did you take xnto account the rotary inertia of the rather large top mass? 

V The direct solution of the beam governing equation would probably give more 
accurate results. The main emphasis of this paper ^ however, has been the 
^i^pl^ methods , such as the Rayleigh^Eitz method as a hand 
calculator code, such as program TUSF, 

2) For the two towers presented in the paper, a Rohr 55V and a Rohr 25G, the 

rotary inertia effects were neglected in the analysis results. These effects 
could be^ easily taken into account, in the simple methods, by adding the 
appropriate term to the kinetic energy, 

\ ^«W. Perkins 

How do you determine a priori the location of nodes for higher mode analysis using 
Rayleigh’s method? 

complex mass and stiffness distributions, the precise mode shape, 

4 priori, is difficult if not impossible to determine , The nodal points are there- 
fore unknown, Jn the Rayleigh or Rayleigh Ritz procedure we choose a shape which 

true mode shape. If the resulting frequency estimate, 
bhis assumed shape, is inaccurate compared to test results, a better mode 
shape approximation is needed. 
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